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Abstract: Not yet. ALEPH's recent exclusion limit employs an aggres- 
sive determination of theoretical uncertainties using a simplified application 
of the Bayesian method. The validity of their analysis can be evaluated 
by its further implications, such as contradicting the existence a b quark 
and requiring relations between hadronic event-shape observables which are 
not observed. Traditional error estimation methods result in a much larger 
estimate for the theoretical uncertainties. This puts the ALEPH and also 
Csikor-Fodor limits at the ~ 1 a level for the very light gluino scenario. A 
recent astrophysical result implies direct searches will be more difficult than 
previously anticipated, adding to the importance of reducing the QCD un- 
certainty in predictions sensitive to indirect effects of light gluinos. Some 
possible indications in favor of a light gluino are noted. 



^Invited talk at Rencontres de la Vallee d'Aoste, La Thuile, March 1997. Research 
supported in part by NSF-PHY-94-23002. 



All the LEP experiments have addressed the question of constraining 
QCD color factors by high-statistics data on hadronic event shapes in Z° 
decay. See, e.g., |1|, 0, ||, Q and references therein. The agreement with QCD 
is entirely adequate, but the sensitivity to the effective number of flavors 
has not been sufficient up to now to rule out the presence of a light gluino 
which would increase the effective number of flavors from ~ 5 to ~ 8. The 
strategy of the new ALEPH analysis is to fit both 4-jet angular distribu- 
tions and the differential 2-jet rate^to simultaneously constrain as{mz) and 
Hf, the effective number of flavors. Since every hadronic event contributes 
to the differential two jet rate, it is well-determined statistically. It is a 
first-order observable, i.e., its leading contribution is 0(0;^), so it is sensi- 
tive to the running of a^. Thus it feels nj through the coefficient of the f3 
function, bo = — ^Uf^. The various 4-jet angular distributions are 

second-order variables, i.e., at the parton level their leading contribution is 
0(0;^). Particular angular correlations enhance the sensitivity to 4-fermion 
final states in comparison to the dominant qqgg final states. However 4-jet 
angular distributions are less precisely determined statistically than is the 
differential 2-jet rate. The nice observation of ALEPH is that the shapes 
of the error ellipses from the two types of measurements tend to be com- 
plementary]^, so combining the constraints could produce a more precise 
determination than either one alone. The final result quoted by ALEPH 
from this analysis is0]: as{mz) = 0.1162 ± 0.0012(stat) ± 0.0040(s?/st) and 
Uf = 4.24 ± 0.29{stat) ± 1.15{syst). The crucial question is whether their 
estimate of the systematic uncertainty is realistic or overly optimistic. 

The fundamental issue is ALEPH's approach to assigning quantitative 
systematic error bars. For each source of systematic uncertainty considered, 

^For each event, the value of the jet-definition parameter ycut at which the event changes 
from being a 2-jet to a 3-jet event is determined. This quantity is called 2/3. The differential 
2-jet rate {D2) is We discuss ALEPH's fit B, which fixes to the QCD value. 

■^For a more detailed discussion of this point see 
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several variations on the model are fit to the data, resulting in a determi- 
nation of as{mz) and Uf, with some for the fit. They adopt what they 
characterize as a "Bayesian point of view", which they describe as follows]^]: 
"The Bayesian idea is that a priori all models can be considered equally 
well suited for usage in the analysis, but from a bad it is deduced that 
the a posteriori probability of such a model is low, and therefore this model 
should get a small weight when estimating the actual systematic error. ...the 
systematic error corresponds to the increase in by one." If all possible 
models, or at least a complete unbiased sample of models is explored by the 
analysis, the above ansatz should be valid on average over many experimen- 
tal fits. The question which must be addressed is whether the ansatz remains 
valid when restricted to the class of model variations considered by ALEPH0. 

It is intuitively plausible that when data are being modeled from first 
principles with a correct theory, and only the parameters of the theory re- 
quire specification, that the prescription outlined above would give a correct 
estimation of errors. Moving away from this ideal but rarely realized case, 
one may have a situation in which the form of some model function is known 
theoretically (e.g., the power-dependence of pseudoscalar meson masses on 
the light-quark mass in the rrig —>■ limit of a lattice gauge theory calcula- 
tion) and only some parameters of the function are unknown. Here it is also 
possible to rely on the Bayesian analysis as described above, for aspects of 
the modeling which refer to the known functional form. 

However in the situation most commonly encountered in practice, the 
form of functions appearing in the model are not known from first princi- 
ples. Examples are the functional forms which are used to describe parton 
distribution functions or parton hadronization. When, as is usually the case, 
functional forms are chosen for simplicity and convenience, and parameters 

^For a discussion of Bayesian methodology, see ref. [|| . Note that the formulae used in 
ref. 0] to define la systematic errors were developed by ALEPH and a derivation is not 
available in the literature. G. Dissertori, private communication. 
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tuned to fit a diverse collection of data, naive application of the Bayesian 
approach can lead to a completely wrong estimate of the actual theoretical 
uncertainty. A recent case in point is the production of high jets ob- 
served by CDF which was anomalously large in comparison to predictions 
of QCD with the then-standard gluon distribution functions0. With the 
standard functional form of the gluon distribution function, the QCD pre- 
diction seemed very strongly constrained by other data. However only a mild 
generalization of the form was needed to fit all data including the high-p^ 
jetsi. 

The most difficult situation of all is when one must model both the con- 
ventional theory and the theory with some new degree of freedom. In this 
case, to completely survey the model space as required for a correct analysis, 
requires modeling and fitting all the physical observables used to constrain 
the models, with and without the new degree of freedom. Here, since the 
form and parameters of hadronization models are tuned to agree with an 
enormous body of data, the correct application of Bayesian principles mini- 
mally requires generalizing Herwig and JETSET to include light gluinos and 
varying all the parameters of the hadronization models to fit all data. To 
the extent that the underlying physics (e.g., non-perturbative QCD) is well 
or poorly understood, the generalization of the modeling to new degrees of 
freedom will or will not provide an adequate framework for the Bayesian 
analysis. 

In addition to these non-perturbative issues, the perturbative QCD pre- 
dictions at the parton level depend on the renormalization scale /i, the match- 
ing scheme used to obtain resummed predictions, and also on the parameter 
ijcut and scheme used to define 4-jet events. The problem of fully exploring 
the model dependence associated with these parameters is discussed below. 

Given that the full implementation of the Bayesian approach is so de- 
manding, we can ask whether it might still be valid in practice when applied 
in the simplified way adopted by ALEPH. Independent experimental evidence 
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will allow us to answer this question. 

Three sources of uncertainty dominate the ALEPH estimate of their sys- 
tematic error: 

• The uncertainty in the theoretical prediction for D21 the differential 
2-jet rate, due to sensitivity to renormalization scale, /x. 

• The uncertainty in the theoretical prediction for D2 due to sensitivity 
to uncontrolled aspects of the resummation of large logarithms of y-^. 

• The modeling of the hadronization of the partons. 

Two other uncertainties in the theoretical predictions which may be impor- 
tant but which were not considered are: 

• The uncertainty in the theoretical prediction for the 4-jet angular dis- 
tributions due to sensitivity to renormalization scale, 11. 

• The uncertainty in the theoretical prediction for the 4-jet angular dis- 
tribution due to large logarithms of ?/cnt, i-e., the dependence of the 
theoretical prediction on jet definition. 

Two schemes for treating the effects of subleading logarithms of in the 
prediction for the differential 2-jet rate were considered, R matching and log- 
R matching. If subleading logs and higher order corrections are unimportant, 
the schemes should give identical results and should exhibit only weak scale 
dependence. Fig. || from |Q shows their fit results for the two schemes, along 
with the for each fit, as a function of Inf = In^^'^/m?^)- It can be seen 
that the extracted value of n/ is sensitive to /i, in both schemes. Since Uf is 
a physical quantity, it cannot depend on /i or scheme. 

If one knew that the true prediction for the functional dependence on 1/3 
of the differential 2-jet distribution were given by the functional dependence 
predicted by one of the schemes, for some particular value of /i, then one could 
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obtain the uncertainty associated with the /x dependence by the Bayesian 
prescription, and would get the correct mean value of the physical parameters 
by using the scheme and fi value with the smallest x^- ALEPH assumes this 
is the case and takes the 1-sigma "theoretical error" (scale uncertainty) to 
be^ (c.f., formulae (15,16) of ref . i) 

A^^nj = C max[,] i ~f ' ==) ; C = max (l, yx^/Aw) . (1) 

y y'max(l, Ix* — ^ ^ 

i labels the fit parameters of various models, with * designating the best fit. 
Table |I|, taken from Table 3 of ref. , shows for the two schemes the best-fit 
and = ±1 fits. Rather than treating i as a continuous variable regarding 
the /X dependence shown in Fig. ^ they only consider the discrete choices 
listed in Table |; applying eqn (P, they find A*'' = l.Ol.Q 

The best fit (log R matching with Inf = —1.3) gives Uf = 3.68. However 
ALEPH averages the best-fit results of the log R and R matching schemes 
and quotes nj = 4.24 ± lM{th) ± 0A5{had) ± 0.27(rfe) ± 0.29{stat). 



Table 1: Scheme and scale variations in ALEPH Fit B analysis, N^of = 73. 



Variations 


as(Mz) 




x' 




Theoretical Prediction 


nominal: logR, In / = —1.3 


0.1154 


3.68 


78.5 


0.52 


logR, In / = -1.8 


0.1182 


4.08 


79.8 


0.41 


logR, In / = -0.6 


0.1130 


3.09 


79.6 


0.74 


R,\n f = -0.5 


0.1210 


5.81 


83.3 


0.78 


i?. In / = 0.0 


0.1175 


4.88 


81.6 


1.00 


R,\n f = 0.7 


0.1141 


3.57 


83.0 


1.42 



The traditional prescription for estimating the ±lo" systematic error from 

scale uncertainty is to take it to be the range of Uf obtained by varying 

/x/\/i over the range [1/2,2]^, i.e., Inf between ±1.4. From Fig. ^ this 

^Applying eqn (|l|) but considering i to run over all /i gives A*'' = 1.1. 
^In the recent analysis the range [1/3,3] is also considered. 
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gives Hf = 4.88 ± 2.5{th), for the R-matching scheme, and Uf = 2.7 ± l{th) 
for the logR scheme. The conservative conclusion is that the dependence of 
measured distributions on is sensitive to truncation of the perturbation 
series (reflected in the sensitivity to /i) and also to subleading logarithms of 
Us (reflected in the sensitivity to matching scheme) and therefore one cannot 
draw strong conclusions until the theoretical predictions are improved. 

Fortunately, we have two pieces of independent empirical evidence to test 
the validity of ALEPH's more aggressive estimate of the scale and resumma- 
tion scheme uncertainty. 

First of all, consider the procedure for fixing the scale uncertainty. As 
noted above, the underlying assumption is that there is some value of fi for 
which the net effect of neglected higher orders and subleading logs is minimal 
so that the functional dependence of the predicted distribution on 7/3 gives a 
good representation of the complete result. This amounts to the principle of 
experimental optimization: fixing /i to the value fi = /igo which minimizes 

for that observable. Unless the differential two-jet rate is a lucky accident, 
the same principle should apply to any event-shape distribution (e.g., thrust, 
energy-energy correlation, etc). If using /iEo does in fact subsume the effects 
of neglected terms, the values of as{mz) obtained from different event shape 
distributions, each at its own fi^o, should be the same up to statistical and 
systematic errors other than those associated with scale uncertainty. Burrows 
has analyzed various scale-fixing proposals ||1CI|| and finds that the dispersion oi 
values of as{mz) obtained with the experimentally-optimized-scale prescrip- 
tion is no better than the dispersion from other scale-fixing prescriptions, as 
evident in Fig. || from ref. reproduced here. That study used 0{a1) 
rather than resummed predictions. However from Figs 34-37 of ref. [|ll[] one 
can see that the dispersion in as{rnz) determinations using ^ieo is large in 
this well, even restricting to matching schemes and observables for 

which l^dof < 1- Thus it is not true that the effects of neglected higher 
order corrections are effectively included by choosing the /i which minimizes 
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Next consider the problem of resummation scheme ambiguity. The log-R 
matching scheme gives the best-fit, with a of 78.5 compared to 81.6 from 
the R-matching scheme. Therefore, if the Bayesian procedure were being 
consistently implemented, the log-R matching scheme would be a posteriori 
identified as the model which gives the best estimate of the physical pa- 
rameters. Within the set of models considered by ALEPH, their analysis 
procedure thus gives Uf = 3.68 ± 0.29(stat) ± 1.15{syst) as discussed above. 
The ±0.29 statistical error means that repeating the experiment many times 
will yield a value of rij in the range 2.93-4.42 in 99% of the trials. The log-R 
matching scheme prediction simply does not tolerate rif > 4.5, for any fi, as 
apparent from Fig. |l|. Thus strictly following Bayesian reasoning requires 
reject the log R matching scheme for the D2 observable because it contra- 
dicts our a priori knowledge that > 5. This phenomenon illustrates that 
Bayesian reasoning can only be applied safely when all relevant experimental 
information is included in the fit and a sufficiently complete set of theoretical 
models is being fit. Otherwise, one can wind up in a local rather than global 
minimum in the space of possible models. 

Hadronization errors are similarly not amenable to the Bayesian method, 
unless much more general hadronization models are used. In particular, the 
parton shower modeling must include light gluinos and the hadrons contain- 
ing them, if one is to consistently determine whether a model with light 
gluinos can fit the data. In the absence of such an undertaking, one can ap- 
ply the traditional method of trying various hadronization models and seeing 
if the resultant values of n/ change. OPAL|jl| investigated the effect on 4-jet 
angular distributions of changing parameters in a given hadronization model 
by 0(10%) and found it produced shifts in Uf of order 1 unit (see Table 1 
of ALEPH did not try changing parameters in a given hadronization 
model, but reports that using Herwig rather than JETSET to model the dif- 
ferential 2-jet rate leads to a best-fit value of rij = 6.21 with only a slight loss 
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of fit quality {x^/Ndof increases from 81.6/73 = 1.12 to 91.6/73 = 1.25). The 
fact that such a large shift mnj accompanies a change of model should have 
been followed up by varying and changing matching scheme, to discover 
the best fit which can be obtained with the new hadronization model. 

A final area where there may be additional theoretical uncertainty which 
cannot yet be quantified is the question of the sensitivity of the 4-jet distri- 
butions to the scale and to details of the jet definition {ycut and scheme). 
The 1-loop correction to 4-jet matrix elements was not available when the 



ALEPH analysis was started. However it has recently been determined [|12 



and the correction to the angular dependence is small. This can be seen in 
Fig. 1^ from ref. |jl3[ which shows the Bengtsson-Zerwas angular distribution 
at tree level and 1-loop for nj = 5 and at 1-loop for nj = 8, compared to 
ALEPH data. Since 1-loop corrections to the angular distributions are small, 
II dependence of the angular distribution is evidently not a problem. At the 
same time. Fig. |^ illustrates that light gluinos make only a very small change 
in the angular distribution, which is why the nj extraction is so sensitive to 
corrections to the perturbative prediction. 

No resummed calculation exists for the 4-jet angular distributions, but 
the effects of resummation may be large. First of all, L3 found consider- 
able sensitivity to t/cut in their color factor analysis P]. Secondly, differential 
distributions such as D2 which depend on defining jets, for which the ycut 
dependence has been determined, display a significant sensitivity to i/cut- Fi- 
nally, the sensitivity of n/ to hadronization model shows that low pt aspects of 
the distribution of hadrons relative to the jets affect the angular distribution. 
Thus soft and colinear gluon radiation may modify the angular dependence 
found in a fixed order calculation. 

It is desirable to estimate what result ALEPH would have reported had 
traditional error estimation methods of earlier experiments [l|, ^ ^ been used. 
We can only give a lower bound on the systematic error, since ALEPH did 
not check the sensitivity of its nj extraction to variation of the parameters 
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of the hadronization model or choice of i/cut for the 4-jet angular distribu- 
tions, while the former made a significant contribution to, e.g., the OPAL 
error We must use R matching, since the log R matching scheme is ex- 
cluded because it contradicts the existence of the b quark. We take cr*^ = 2.5 
because this reproduces the nf range obtained for 1/2 < fi/Mz < 2. Com- 
bining quadratically with the ±0A5{had) ±0.27(c?e) ±0.29{stat) obtained by 
ALEPH, gives Uf = 4.88 ± 0.29{stat) ± 2.57{syst). Using the PDG prescrip- 
tion for obtaining limits when part of the experimental range is unphysical 
{uf < 5), this implies the 95% cl limit Uf < 9.95, or the 77% cl limit nj < 8. 

It is interesting that the ±2.57{syst) obtained from the ALEPH combined 
fit to D2 and the 4-jet angular distributions by traditional uncertainty esti- 
mation procedures is actually a larger than obtained using the 4-jet angular 
distributions alone|[|, 0, |[. The reason for this is that D2 is actually very 
sensitive the /i, even after resummation, while the 4-jet angular distribution 
is not, so that the ALEPH combined analysis only results in a more precise 
determination of Uf if the fi range can be limited. However the results of 
ref. |10| show that ALEPH's attempt to limit the fi range by embracing the 
"experimental optimization" prescription is not valid. 

Csikor and Fodor (CP) receintly studied the effect of light gluinos on the 
running of By restricting themselves to extracted via 3-loop 

perturbative predictions for Rhad in t decay and in e'^e" collisions at the 

and below, they largely avoid the hadronization and scale uncertainties 
of the hadronic event shape analyses. The fit is mainly constrained by the 
values as{mr) = 0.335 ± 0.053 and as{mz) = 0.123 ± 0.006 because the 
data at intermediate values of Q are statistically much weaker. However 
the error on ^^(mT-) has been argued[0 to be of order twice as large as 
used by CP. Thus in their final analysis for the very light gluino case, CP 
do not use asifrir) at all. They find a very light gluino is disfavored only 
at about 70% cl. The 99.97% cl limit they quote after combining with the 
ALEPH limit mainly refiects the ALEPH number, which we have argued 
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underestimates the theoretical uncertainty and does not correctly deal with 
the manifestly unphysical log R matching. Modifying the ALEPH result to 
Uf = 4.88 ± 0.29{stat) ± 2.57{syst), as obtained using standard methods of 
error estimation, leads to the conclusion that very light gluinos cannot be 
excluded with either the ALEPH or Csikor-Fodor methods until QCD uncer- 
tainties, which are at this point mainly non-perturbative, can be reduced. 

Before closing, I remark that SUSY model-building developments and re- 
cent cosmological constraints have increased the importance of such indirect 
constraints on light gluinos. Several gauge-mediated SUSY-breaking models 
have recently appeared in which the gluino is lighter than the other gauginos 



and i?-parity is conserved ||16||. In such scenarios the only decay mechanism for 
the lightest gluino- containing hadron (the glueballino, denoted R^) is through 
gravitino emission, if kinematically allowed. Thus the may be stable, or 
so long-lived that direct detection through its decay may be impossible. Fur- 
thermore, direct detection of decaying i?°'s may also be more difficult than 
anticipated in models in which the photino as well as the gluino is ap- 
proximately massless at tree level. In this case, radiative corrections give 
the gluino and photino masses of ~ 0.1 GeV and < 1 GeV respectively, and 
the i?" mass should be 1.3-2.2 GeV|T^. A new analysis|TB| of the photino 
relic abundance shows that 771(7) ^ m{R^) /I.8. Combining this with the 
1.3-2.2 GeV mass range of the R^ leads to discouraging prospects for direct 
searches for the decay i?" —>■ 7r+7r~7 over the full r = m{R^)/m{'j) range 



in which photinos provide the observed dark matter: 1.4 < r < 1.7|]T8|. The 
point is that the invariant mass of the dipion system cannot be larger than 
m{R^) — m(7) = m{R^){r — l)/r, while to remove the vr^vr" and other 

background one wishes to require m{iT~^TT') convincingly above rriK- Also, 
the low Q value makes it more difficult to use the pt of the tt+tt^ pair as a 
discriminant against background. 

Thus even if there are light unstable gluino-containing hadrons, they may 
be very difficult to observe through their decays. This increases the impor- 
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tance of reducing the (mainly non-perturbative) QCD errors in theoretical 
predictions for as running and for e^e~ event shape observables. When lat- 
tice gauge theory can predict the bottomonium spectrum without the use of 
non-relativistic or quenched approximations, the nf required to obtain pre- 
cision agreement for the level splitting may provide the cleanest indirect test 



of all. Direct searches for quasi-stable R-hadrons as discussed in ||19| , |20| , [17 
should be considered. 

My La Thuile talk ended with a brief discussion of three tantalizing phe- 
nomena which could be naturally explained by light gluinos but are seemingly 
difficult to explain otherwise. Here, I merely list them in order of decreasing 
experimental robustness and leave the reader to consult the references given 
for details and further references. 

• The existence of the isosinglet pseudoscalar ry(1410) whose properties 
match those expected for a gluinoball or glueball and whose mass agrees 
with that predicted for a gluinoball but is much lighter than the lattice 



QCD and sum rule predictions for a glueball |21 



Ultra-high-energy cosmic rays which violate the Greisen-Zatsepin-Kuzmin 



bound 22, M 



• The peak observed in dijet mass pairs by ALEPH in e^e annihilation 
at 133, 161 and 172 GeV. 
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Figure 1: Fit results with the R (sohd) and logR (dashed) matching schemes 
function of Inf = Infi'^/s, from @]. 
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Figure 2: (a) values of as{mz) from QCD fits to various event shape dis- 
tributions using "physical" (/i = Mz) - solid circles - and experimentally- 
optimized scales - open circles. For each observable the shaded region indi- 
cates the total uncertainty estimated in Ref. 
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dominated by the contri- 
bution from wide variation of the renormalization scale, (b) shows similar 
comparisons for other scale-setting prescriptions: Principle of Minimal Sen- 
sitivity (solid circles), Fastest Appa^^^t Convergence (solid triangles) and 
Brodsky-Lepage-Mackenzie (open squares). Taken from ref. [p!0| . 




Figure 3: Bengtsson-Zerwas angular distribution at tree and 1-loop level for 
ri/ = 5 (dotted and solid) and at 1-loop for n/ = 8, compared to ALEPH 
data, from |jl3[. 
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